A multi-dimensional numerical model has been developed for the unsteady state oscillatory combustion of solid propellants subject to acoustic pressure disturbances. Including the gas phase unsteady effects, the assumption of uniform pressure across the flame zone, which has been conventionally used, is relaxed so that a higher frequency response in the long flame of a double-base propellant can be calculated. The formulation is based on a premixed, laminar flame with a one-step overall chemical reaction and the Arrhenius law of decomposition with no condensed phase reaction. In a given geometry, the Galerkin finite element solution shows the strong resonance and damping effect at the lower frequencies, similar to the result of Denison and Baum. Extended studies deal with the higher frequency region where the pressure varies in the flame thickness. The nonlinear system behavior is investigated by carrying out the second order expansion in wave amplitude when the acoustic pressure oscillations are finite in amplitude. Offset in the burning rate shows a negative sign in the whole frequency region considered, and it verifies the experimental results of Price. Finally, the velocity coupling in the two-dimensional model is discussed.